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First observations of XcO, Xa, and Xc2 decays to tt^tv^ K''^Kg and K'^ K~ KgKg, as well as ip{2S) 
decay to tt^tv' KgKg are presented. The branching fractions of these decay channels are determined 
using 14 X 10*^ V'(2S') events collected at BESII/BEPC. The branching fractions of XcO, Xc2 ^ K°gK°g 
are measured with improved statistical precision. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 



§ The denote charged pions or kaons. 
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I. INTRODUCTION 

Experimental data on charmonia and their decay properties are essential input to test QCD models and QCD 
based calculations. The importance of the Color Octet Mechanism(COM) in radiative decays of the T 0, J/V" 
production in inclusive B decays 0, as well as inclusive decays of P-wave charmonia has been emphasized for 
many years. Recently, QCD predictions of two-body exclusive decays of P-wave charmonium with the inclusion of 
the COM have been made 'F, v\ and compared to previous measurements 0, Q • More experimental data of two- and 
four-body exclusive decays of P-wave charmonia with improved precision are important for further testing this new 
QCD approach including the effect of the COM. 

In this paper, results on ipi'^S) and Xcj {J = 0, 1, 2) two- and four-body hadronic decays with inclusion of a pair of 
K'g mesons are presented. This analysis is based on 14 x 10^ i'i'^S) decays collected with BESII at the BEPC e+e^ 
Collider. A sample of 6.42 pb^^ data taken at 3.65 GeV is used for continuum background studies. 

II. BES DETECTOR 

The BESII detector is described elsewhere Charged particle momenta are determined with a resolution of 
fp/p = 1.78%\/l -I- (p in GeV/c) in a 40-layer main drift chamber (MDC). Particle identification is accomplished 
using specific ionization (dE/dx) information in the drift chamber and time-of-flight (TOF) information in a barrel- 
like array of 48 scintillation counters. The dE/dx resolution is cr^E/dx = 8%; the TOF resolution is (Ttof=200 ps 
for hadrons. A 12-radiation-length barrel shower counter (BSC) measures energies of photons with a resolution of 
as/E = 21%/ Ve {E in GeV). 

III. MONTE CARLO SIMULATION 

A GeantS based Monte Carlo, SIMBES which simulates the detector response, including interactions of sec- 
ondary particles in the detector material, is used to determine detection efficiencies and mass resolutions, as well as to 
optimize selection criteria and estimate backgrounds. Under the assumption of a pure El transition, the distribution 
of polar angle 6 of the photon in ip{2S) JXcJ decays is given by 1-t-fccos^ 6 ITl| with fc = 1, — ^, and for J = 0, 1, 
and 2, respectively. The angular distributions for Kg mesons from XcO.2 — KgKg decays are produced according 
to the model of Xcj MM where M stands for a 0~ meson. Angular distributions for daughters from other 
decays are generated isotropically in the center-of-mass system of the iJj{2S) or Xcj- 

IV. DATA ANALYSIS 

To be regarded as a good photon, a shower cluster in the BSC must have an energy deposit of more than 50 
MeV and at least one hit in the first six layers of the BSC. To remove soft photons emitted by charged particles, the 
differences of azimuthal angles, d0, and z coordinates at the first layer of the BSC, dz, between good photons and each 
charged track must satisfy either a loose requirement (selection-A: dcj) > 10° or dz > 0.3 m) or a tight requirement 
(selection-B: dcj) > 20° or dz > 1.0 m). Here the z coordinate is defined to point in the positron direction. 

Each charged track is required to have a good helix fit. For final states containing charged kaons, particle identifi- 
cation is required; usable particle identification information in one or both of the MDC (dE/dx) and TOF subsystems 
is necessary. A particle identification x^ is calculated for each track for the pion, kaon or proton hypotheses using 
this information, and the associated probability prob is determined. A track is identified as a kaon, if the probability 
of the track being a kaon prob(K) > 0.01; otherwise it is regarded as a pion. For final states containing only pions, 
no particle identification is done and all tracks are assumed to be pions. 

Each event is required to contain two Kg mesons. The reconstruction of the decay Kg tt+tt^ and related checks 
are described in detail elsewhere A Kg candidate must satisfy \M^+^- — Mj^o\ < 20 MeV and have a decay 

length transverse to the beam axis Rxy > 0.3 cm. The Kg sideband sample, used for background estimation, is 
selected with one tt'^'tt" pair within the Kg mass window and the other pair in the Kg mass sideband region defined 
by 40 MeV < \M^+^- - Af^gl < 60 MeV. 

Four constraint (4C) kinematic fits are performed on the selected events for the following decay modes : (1) ip{2S) 
-fK°gK°g, (2) i}{2S) jTT+Tr~K'^K°, and (3) ij{2S) -fK+K-K°gK'^. The fits are made to each combination of a 
good photon and two Kg candidates in an event, the combination with the minimum xlc is selected, and the xlc is 
required to be less than 35. The associated probability prob^c is calculated. 
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Background from Tr~^TT~J/ijj decay is removed by calculating the mass recoiling, MrecoU, against all pairs 

of oppositely charged tracks, assuming them to be pions, and requiring \Mrecoii — Mj/^\ > 25 MeV. Background 
contamination from continuum production is found to be negligible for all decay channels. 

An unbinned maximum likelihood method is used in fitting the signal for all decay channels except 4^(23) — *■ 
h^h~ KgKg. The branching fractions of ip{2S) IXcj [J — 0,1,2) needed in the measurement are taken from 
Particle Data Group (PDG) tables |3]. 



t/.(2S) ^ -iK%K% 



The decay ip{2S) — > jK^Kg has one photon plus a pair of Kg candidates. The event should have four charged 
tracks with total charge zero. The loose photon selection, selection-A, is applied because of the low background in the 
channel. The KgKg invariant mass distribution of the selected events is shown in Fig. ^ A few Kg sideband events 
survive the selection, which is consistent with the low background observed in Fig.n(a). No background is expected 
from tp{2S) -jXcj with Xcj 2(7r+7r") for J = 0, 1, 2 and ip{2S) jxa with xa 
analysis of simulated MC events. 

The KgKg invariant mass distribution is fitted with two Breit-Wigner resonances for Xco and Xc2, each convoluted 
with Gaussian resolution functions, plus a second order polynomial background. The Xc0,2 widths in the fitting are 
fixed to their PDG values 3] • The resulting fit is shown in Fig. ^ (b) . Including the Xci resonance in the fit yields 
zero events for the CP violating decay Xci ^s-^s- 



KgK^TT^ according to the 
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FIG. 1: Distribution of KgKg invariant mass of ip{2S) 'jKgKg candidates, (a) Points with error bars are data, and the 
liistogram is sideband background, (b) Points with error bars are data, and the solid line is the fit described in the text. 
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The ipi^S) — * 'yn^TT^ KgKg decay channel contains one photon and six charged tracks with total charge zero. The 
requirements here are similar to the previous case, but there are two additional pions. Background from ir/K misiden- 
tification is suppressed by the requirement prob4c{'fn'~^TT~ K^Kg) > prob4c{jK~^K~KgKg). The tt^tt~ K^Kg in- 
variant mass distribution for selected events is shown in Fig. [5] 

In Fig. [3 there are two kinds of background in the mass region between 3.0 and 3.64 GeV/c^: (1) background 
corresponding to Kg sidebands, and (2) iIj{2S) decays and XcJ decays different from the signal channel, where the 
decays also include a pair of Kg mesons. Studies with Kg sideband events for both data and MC show that Kg 
sideband background from wrong combinations of tt+tt" is slightly enhanced in the Xcj signal region. MC studies 
show that the smooth background spread over the whole mass region from (2) results mainly from the following decay 
channels: (a) ^(25') ~> -fXcj with Xd 3(7r+7r-) and Xd K+K-K^K^, (b) ^(25) ^ t:^t:+i:^ K^K^, and 
(c) ipi^S) — > ujKgKg with Lu 7r+7r~7r°. Background events in the high mass region above 3.64 GeV/c^ in Fig. [21 
are from V'(2>S') Tr~^Tr~ KgKg decays combined with an unassociated low energy photon. 

The TT'^TT~ KgKg invariant mass distribution between 3.0 to 3.64 GeV/c^ is fitted with three Breit-Wigner resonances 
XcJ {J — 0,1,2), convoluted with Gaussian resolution functions, plus a second order polynomial background. The 
widths of the XcO,i,2 resonances in the fit are fixed to their PDG values. The fit is shown in Fig. 12 The numbers 
of events in the three peaks determined from the fit include signal and Kg sideband background, which is somewhat 
enhanced in the regions of the peaks. The Kg sideband sample for data is fitted with a fake signal shape, found 
by fitting the MG Kg sideband sample, plus a second order polynomial background. The numbers of sideband 
background events, 5.3, 0.6 and 5.5 for Xco^ Xci and Xc2 respectively, are then subtracted from the total numbers of 
events in three peaks. 




l\/l(nnK°sK°s) (GeV/c^) 



FIG. 2: Distribution of ■k'^-k~ KgKg invariant mass for ip{2S) — > jtt^tv' KgK^ candidates. Points with error bars are data. 
The light shaded area in (a) is background simulation, where some unknown branching ratios are normalized to agree with the 
overall XcJ background level, and the dark shaded area is Kg sideband. The solid line in (b) is the fit. 
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C. ip[2S) jK+K-K°sK°s 

The tp{2S) 'fK^K^KgKg decay has the same topology as ip{2S) 'fn^TT^ KgKg, and thus it is subject to 
similar event selection criteria except for the kaon identification requirement for two of the charged tracks. First, the 
KgKg pair is searched for under the assumption that all charged tracks are pions. Kaon identification is only done for 
the two charged tracks remaining after reconstruction of the KgKg pair. We also require prob^ci'jK'^K^KgKg) > 
prob4c{l'^^''^~ KgKg) for the 4C kinematic fit probabilities to suppress contamination from ^{2S) 'jTr^n^ KgKg 
decays. The K'^K~KgKg invariant mass distribution for selected events is shown in Fig. 13 

As seen from Fig. O only one event survives from the Kg sideband sample for data. MC events for the following 
possible background channels are generated: (1) ■ip{2S) JXcJ with Xcj 3(7r+7r~) and tt'^tt" K^Kg, (2) tp{2S) 
7T~^Tr~ KgKg, and (3) ip{2S) wKgKg with uj tt+tt^tt'^. However, no event from these background channels 
survives the selection criteria. Another study with a large sample of simulated ip{2S) — > anything 0| shows that 
negligible background comes from decays of iIj{2S) — > (t>K*^K^ tt^ K~ K^K^. 

The K^K^K^K^ invariant mass distribution is fitted with three Breit-Wigner resonances, XcJ {J — 0,1,2), 
convoluted with Gaussian resolution functions, plus a flat background. Because of low statistics in the signal region, 
not only the widths and mass resolutions for the Xcj ("^ = 0, 1, 2), but also the masses of the Xci and Xc2 in the fitting 
are fixed to their PDG values. The fitting results are shown in the Fig. 13 




3 3.1 3.2 3.3 3.4 3.5 3.6 3.7 



M(lCKrK°sK°s) (GeV/c^) 

FIG. 3: Distribution of K+R-R^K^s invariant mass of ?/>(25) -jR+R-R^sRs candidates. Points with error bars are data, 
and the histogram is sideband background. The solid line is the fit. 



D. V(2S) ^ h+h-R'^R'^ 

The selection of ip{2S) — > h'^h~ KgKg decays requires six charged tracks with total charge zero and no good 
photon in the event, as defined above. Good photons are rejected with the tight selection, selection-B, in order to 
gain higher detection efficiency for signal events. The Kg reconstruction uses all combinations of oppositely charged 
tracks assuming all tracks are pions. To further suppress background of ip{2S) radiative decays, a requirement on the 
missing momentum of six charged tracks is employed: Pmiss < 80 MeV. The two charged tracks and h~ recoiling 
against the Kg pair are assumed to have the same mass m. Using energy-momentum conservation, the mass squared 
is calculated from 

where E = M^(2S) ~ Ef^oj^o, and Pfi± is the momentum of /i+ or . The distribution of for selected events is 
shown in Fig. 01 The peak at low mass is consistent with tt+tt^; there is no evidence for K^K^. 

Two events from the continuum data sample survive the above selection and their effect will be included in the 
systematic error. No background is found in MC studies of the following decay channels : (1) '0(25') IXcJ 
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with XcJ 3(7r+7r-), Tr+Tr-ifO ^ , and K+R-RlKl and (2) V(2S') ^ ^jK^KI with w ^ tt+tt-ttO. Background 
estimated using the sideband data is subtracted from the observed number of signal events. A MC study shows 
that the shape of the charged pion signal in the spectrum is well described by a Gaussian function, and its mean 
and resolution are consistent with data. The spectrum is fitted with a Gaussian signal function and a flat background 
using a binned maximum likelihood fit where the resolution is fixed to the MC determined value. The fitting result 
is shown in the Fig. ^ 
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FIG. 4: Distribution of invariant mass squared of the two remaining charged particles after KgKg selection for ip{2S) — » 
h'^ h~ KgK^. (a) Points with error bars are data. The histogram is the Kg sideband background, (b) Points with error bars 
are the data with the Kg sideband background subtracted. The solid line is the fit. 



E. Systematic Errors 



Systematic errors for the efficiency are caused by differences between data and MC simulation. Our studies have 
determined these errors to be 2% per track for the tracking efficiency, 2% for photon identification, 5% for the 4C 
kinematic fit, and 2.1% for the Kg reconstruction efhciency. A correction factor due to the overestimate of the Kg 
reconstruction efficiency of the MC relative to data is determined to be 95.8%. The change of fitting range and 
background shape function contributes a difference of final results less than 3%. Other systematic errors arise from 
the uncertainties in the total number of ip{2S) events, (14.00 ± 0.56) x 10^ ^^'^ branching fractions for 

Kg 7r+7r~ and ip{2S) — » 7Xcj (J=0,I,2). In ip{2S) tt^tt^ KgKg decay, with two events found in contimmm 
data, an additional error of 7.7% is added. 
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TABLE I: Summary of the fitting results. Errors for the signal yield n^, background ni,, mass M, and mass squared are 
statistical. The detection efficiency e and resolution a for each decay channel from MC are shown. 



Channel 


Us 




(MeV/c^) 


e 

(%) 


(MeV/c^) 




322 ± 20 




3413. 1±1. 2 


7.96 


13.3 


Xcl KgKs 





6.4 ± 2.6 


hxed 


8.50 


12.8 


^ JVO R'O 


DO.i ±0.1 




oooo. ( ±i.o 


0.48 


1 1 o 
11. 


\^ r. ^ ■w + -n-~ 
XcO n TT 


1 ^9 + 1 4 




■3^11 9 + 90 


2.03 


16.8 


Xcl ^ /t /I l\gl\c; 


19.8 ± 7.7 




3,^m 1 ± 6 2 


2.20 


16.4 


Xc2 ^ 7T + 7V- K°sK°s 


57± 11 




3548.2 ± 3.1 


2.04 


17.2 


XcO ^ K+K-K"sK"s 


16.8 ±4.8 




3415.4 ± 6.1 


0.91 


16.1 


Xcl ^ K+K-K%K% 


3.2 ± 2.4 


1.8 ±0.8 


fixed 


1.12 


15.3 


Xc2 ^ K+K-K%K% 


2.3 ± 2.2 


1.8 ±0.8 


fixed 


1.05 


15.9 


Channel 




rib 


m"(10-^) 


e 


a(10-^) 








(GeVVc^) 


(%) (GeVVc*) 




83.2 ± 9.4 




18.0 ±3.1 


2.82 


26.5 



TABLE IL The branching fractions from this measurement, as well as previous results, are listed. The first and second errors 
for the branching fractions BR are statistical and systematic, respectively. 



Channel 


BR{^{2S) ^ 7Xc) X BRixc ^ X) BR[xc ~* X) 

(10-^0 (10-") 


BRpDoiXc ^X)\^ 
(10-*) 


XcO - K'iK^s 
Xcl ^ K%K% 
Xc2 ^ K%K% 


30.2 ±1.9 ±3.3 35.1 ±2.2 ±4.7 
< 0.6 (CL=90%) < 0.8 (CL=90%) 
5.72 ±0.76 ±0.63 8.9 ± 1.2 ± 1.3 


21 ±6 
7.2 ± 2.7 


XcO ^ ^+T^-K^sK"s 
Xcl ^+T^-K%Kl 
Xc2 ^ ^+T^-K%K% 


55.8 ±5.1 ±8.9 65±6±12 
6.7 ±2.6 ±1.1 8.0 ±3.1 ±1.5 
20.7 ±3.9 ±3.3 32.4 ±6.1 ±6.2 




XcO-* K+K-K"sK"s 
Xcl ~*K+K-K%K% 

Xc2^K+K-KlK% 


13.8 ±3.9 ±2.5 16.0 ±4.6 ±3.2 
2.1 ±1.6 ±0.4 2.5 ±1.9 ±0.5 

< 4.2 (CL=90%) < 5.1 (CL=90%) 
1.6 ±1.6 ±0.3 2.6 ±2.4 ±0.5 

< 3.5 (CL=90%) < 5.5 (CL=90%) 




Channel 


BR{iP{2S) X) 
(10-^) 


BRpdg{^{2S) ^ X)\8] 
(10-*) 


iIj{2S) tv^tt'K"sKs 


2.20 ±0.25 ±0.37 





F. Result and Discussion 

Possible resonance structures have been searched for the XcO tt^tt" KgKg final state which is the channel with 
the highest number of observed events. Some excess for inclusive decays of i^*(892)+ Kgir^, /o(1710) KgKg, 
p(770) 7r+7r~ and /o(980) tt+tt" can be seen from the selected events. Insufficient statistics and compli- 
cated structures in these decay modes make it difficult to identify clear signals for two-body decays with interme- 
diate resonances. Efficiencies for final states with resonances, such as 7^*(892)+i^*(892)~, (1430)+ JCn (1430)", 
ii:o*(1430)+i4:2*(143G)-, /o(1370)/o(1710), /o(980)/o(98G), /o(980)/o(22GO) and ifi(127G)''A'° 16] are studied using 
phase space MC events. The averaged difference in efficiency between final states with and without intermediate res- 
onance is estimated to be 7.7%, which is regarded as systematic error in the measurements of the branching fractions 
for the four-body final states. The results of four-body final states h^h~ KgKg in our mesurements include those of 
both non-resonance and intermediate resonance. 

Final results of signal yield and branching fractions for the Xcj(l-P) and ip{2S) two- and four-body hadronic decays 
involving Kg pair production are summarized in TableQ] The masses of the XcJ (J=0,l,2) extracted from the fits are 
also listed. The 90% confidence level (CL) upper limits on the branching fractions in the table are obtained using the 
Feldman-Cousins method [l3|- The branching fractions of Xcj (J=0,l,2) decays to Tr~^TT~ KgKg and K~^K~KgKg, 
as well ipiiS) decay to Tr~^TT~ KgKg are observed for the first time. The branching fractions of XcO and Xc2 decays to 
KgKg are measured with improved precision. 

Decay rates, determined using updated Xcj total widths "s"] and branching fractions for Xcj ~^ tt*^tt*^ , tt+tt" (J = 0, 2) 
and XcJ pp [ J = 1, 2) decays j^], provide support for the COM (see Table ITTTji . According to isospin symmetry, the 
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TABLE III: Comparison of partial widths for Xcj tvtt,KK and pp decays between PDG |2J and the COM predictions. Also 
shown is the result based on this analysis. 



Decay 


V,{PDG) 
in KeV/c^ 


Ti{COM) 
in KeV/c^ 


TT — — 

XcO TT^TT 
Xc2 TT + TT^ 


49.5 ± 6.7 
3.73 ± 0.64 


45.4 f5] 
3.64 f5] 


U U 
XcO ^ TT TT 

Xc2 - 7rV« 


25.3 ± 3.3 
2.3 ± 1.5 


23.5 [5] 
1.93 f5] 


Xcl PP 
Xc2 PP 


0.066 ±0.015 
0.143 ±0.018 


0.05627 [6] 
0.15419 


XcO ^ K+K' 

Xc2 ^ i^+i^- 


61 ± 10 
1.98 ±0.47 


38.6 f5] 
2.89 f5] 


XcO ^ 

Xc2 ^ K^'K'^ 


71 ± 12 (this paper) 
3.76 ± 0.80 (this paper) 





Xcj K^K'^ and K^K decays should have the same partial width. Assuming equal decay widths for Xcj ^s^s 
and K^KI, we find that the partial width of the XcO K^K^ decay estimated using the result obtained in this 
paper is not consistent {2.7a) with the COM prediction for XcO — > K'^K~ , while the agreement between them for 
the corresponding Xc2 decay is within l.lcr. A comparison for the Xcj K^K~ [J — 0,2) decays shows that the 
discrepancy between PDG values and the COM predictions is 2.2a and 1.9(7 for XcO and Xc2 decays, respectively. 

Furthermore, the sum of all known XcO two-body branching fractions is less than 2%. It therefore is important 
to measure more Xcj decay modes, including two-body modes with intermediate resonance and many-body modes, 
because of their large contribution to the hadronic decay width. Theoretical predictions with inclusion of the COM 
for XcJ decays to many-body final states are required for comparison with data. 
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